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Abstract. This paper is devoted to the asymptotics of the density of surfacic states near the spectral edges 
for a discrete surfacic Anderson modeL Two types of spectral edges have to be considered : fluctuating edges 
and stable edges. Each type has its own type of asymptotics. In the case of fluctuating edges, one obtains 
Lifshitz tails the parameters of which are given by the initial operator suitably "reduced" to the surface. For 
stable edges, the surface density of states behaves like the surface density of states of a constant (equal to 
the expectation of the random potential) surface potential. Among the tools used to establish this are the 
asymptotics of the surface density of states for constant surface potentials. 



0. Introduction 

On Z"' {d = di+ d2, di > 0, da > 0), we consider random Hamiltonians of the form 



(N ■ 1 

(N : H^ = --A + V^ 



where 



-A is the free Laplace operator, i.e., — (Au)(n) = X^|„„„|=i u{m); 



■ • Vuj is a random potential concentrated on the sub-lattice Z'^^ x {0} C Z"^ of the form 

O : (0.1) K.(7i,72) = ^ n '7 = (71,72) e x Z'^- = Z'^. 

. if 72 7^ 0. 

o ■ 

^ ' and (<^7i)^^£2di is a family of i.i.d. bounded random variables. For the sake of simplicity, let us 

Q_|. assume that the random variables are uniformly distributed in [a, b] {a < b). 

fH . To keep the exposition as simple as possible in the introduction, we use these quite restrictive assumptions. 
[ We will deal with more general models in the next section. 

■ The operator Hi^ is bounded for almost every uj. It is ergodic with respect to shifts parallel to the 
• . . surface. So we know there exists S the almost sure spectrum of H^^j (see e.g. [14, 23]. 

' For H^^, one defines the integrated density of surface states (the IDSS in the sequel), in the following 
• way (see e.g. [8, 2, 3, 20]): for (p G C^(M), we set 

where IIi is the orthogonal projector on the subspace C^o ^ £'^{1j'^'-^) C £^(Z'^). Here, Sq denotes the vector 
with components ((^ojO^g^di • 

Obviously, equation (0.2) defines the integrated density of surface states only up to a constant. We 
choose this constant so that Ug vanishes below S U Sq where Sq is the spectrum of — ^A. We will see later 
on that, up to addition of a well controlled distribution, is a positive measure. 

One knows that S = (t(— ^A)Usupp(dns) (see [8, 9, 2]. We will study the behavior of Ug at the edges 
of S. To simplify this set as much as possible, we will assume that the support of the random variables 
(u;7j)^^g^di is connected. Under this assumption, we know that 

Lemma 0.1. is a compact interval given by 

(0.3) S = a(-iA,J+ U <r{-^Ad,+ujUl) 

a;oe[a,fe] 

where IIq is the projector on the unit vector 6^ E £^ {Z'^^ ) . 
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This is a consequence of a standard characterization of S in terms of periodic potentials (see [14, 
23]). The assumption that the random variables have connected support can be relaxed; more connected 
components for the support of the random variables will in general give rise to more spectral edges (as in 
the case of bulk randomness, see [16]). For the value of S, two different possibilities occur : 

(1) S = (j(— ^A) + [—a, P] = [—d — a,d + f3] where a = a(o), f3 = /3{b) and a + (3 > 0; this occurs 

• if ^2 < 2 and either a < 0, in which case a{a) > 0, or 6 > 0, in which case /3(6) > 0, 

• if (i2 > 3 and a > or b > bo, where, by (0.3), the thresholds oq and bo are uniquely determined 
by the family of operators (— ^A^j + tnQ)fgiR. 

If a > (resp. /3 > 0), we say that the left (resp. right) edge is a "fluctuation edge" or "fluctuation 
boundary" (see [23]). If a = (resp. f3 = 0), we will speak of a "stable edge" or "stable boundary". 

(2) E = o"(— ^A); this occurs only in ^2 > 3 and if a is not too large, that is, if a G (0, ao]. 
In this case, both spectral edges are stable. 

On the other hand, it is well known (see [24]) that, 

• if ^2 = 1,2, then, for a > 0, (t(— ^A^j — oIIq) = [—^2,^2] U {A(a)}, and the spectrum in [—^2,^2] is 
purely absolutely continuous and A(a) is a simple eigenvalue; 

• if ^2 > 3, there exists > such that 

— if < a < a*^, then, a{—^Ad2 — oHq) = [—^21^2]! and the spectrum is purely absolutely 
continuous; 

— if a = a^, then 

* if d2 = 3,4, then cj(— ^A^^^ ~ ^IIq) = [—d2,d2], the spectrum is purely absolutely contin- 
uous, and —d2 is a resonance for — ^A^fj — oIIq; 

* if (^2 > 5, then cj(— ^A^j ~ '^^^q) = [—d2,d2], the spectrum is purely absolutely continuous 
in [—d2,d2), and — (^2 is a simple eigenvalue for — ^Arfj — oHq] 

— if a > a", then, cj(— ^Aj^g — alio) = [—1^2,^2] U {A(a)}, and the spectrum in [—d2,d2] is purely 
absolutely continuous and A(a) is a simple eigenvalue; 

For the operator — ^A^j + ^Hq, we have a symmetric situation. 

Our aim is to study the density of surface states near the edges of S. In the present case, both edges 
are obviously symmetric. So we will only describe the lower edge. One has to distinguish between the case 
of fluctuation and stable edges. The behavior in the two cases are radically different. 

0.1. The stable edge. As the discussion for lower and upper edge are symmetric, let us assume the lower 
edge is stable and work near that edge. 

In the case of a stable edge, it is convenient to modify the normalization of the IDSS. Therefore, we introduce 
the operator 

Ht = -^A + ti^n^. 

As above, let a be the infimum of the random variables For if G C^(IR), define 

The advantage of this renormalization is that the IDSS 

^^s,norm tliG distributions.! derivcttivG of a positive 

measure. Indeed, for ip E Co^(M), define 

{^,dNs,norm) = -E(tr(ni [^(99) (i^^) - P (if) {H a)]Ui)) 

where 

r+00 

P{^){x) = / ip{t)dt. 

J X 

Clearly, dNs. norm IS independent of the anti-derivative of f chosen to define it; it is a positive measure and 

ah 
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Let n* be the IDSS for Ht- As above, one can define a anti-derivative of n*; denote it by — diV*. Let ^orm 
be the normahzed version of n*, i.e. n^ j^^^^ = n\ — n1. One has 

(0.4) n^^norm + ?^s = <• 

One problem one encounters when studying Ug is that very little is known about its regularity for random 
surfacic models (see nevertheless [21]). Thanks to (0.4), we know that is the difference of two distributions 
each of which is the derivative of a signed measure. So we can take the counting function of dNg as 
dNg = diVs^norm ~ ^Ng is the difference of two measures. Thus, we define its counting function 

(0.5) Ns{E) = r dNs{e). 

J-d 

An obvious consequence of (0.4) is the 
Proposition 0.1. One has 

(0.6) N'^{E) < Ns{E) < N^,{E). 

In section 5.1, we study the asymptotics for N^. As a consequence of this study, we prove 
Theorem 0.1. Assume d2 = 1 or 2. Then, one has 

Fo/(S^i-i) 



NsiE) ~ 

E-*-d 
E>-d 



2VoZ(S'^i~i) (S + (i)i+^i/2 



(^ + d)l+rfi/2 ifd2 = l, 

ifd2 = 2. 



Ui(rfi + 2)(27r)'^i \log{E + d)\ 
where S'^^"^ is the di — 1 dimensional unit sphere. 

If a > 0, this result is an immediate consequence of Proposition 0.1 and of Theorem 1.1 giving the asymptotics 
of the IDSS for constant surface potential (see also section 5.1). If a = 0, one needs to improve upon (0.6) 
as the left hand side of this inequality vanishes making it unusable. This is the purpose of Theorem 1.2. 

When d2 > 3, the situation becomes more complicated and we are only able to use Proposition 0.1 to 
get the two-sided estimate 

^"•^^ ^ (l + al) - s{E + d){E + dy+'^^/^ ■ ^'^^^ - ^ (1 + 6/) 

where C is a positive constant depending only on the dimensions di and d2 (see section 5.1) and 

s{x) = -\x\ 2 , 



\ sup ! 



6»2GT''2 



d - ^ cos(0j) - ^ cos(0^) I dQ2 
i=i i=i 



Here, and in the sequel, the measure dQa [a G {1, 2}) is the Haar measure on the torus T"^", i.e. the Lebesgue 
measure normalized to have total mass equal to one. 

Let us note that, if a < < 6, the inequality (0.7) does not give much information of the actual 
behavior of Ns{E) when d2 > 3. 

0.2. The fluctuation edge. Here, we assume that Eq = ini a{H^) is strictly below —d = inf o"(— ^A). In 
this case, Eq is a fiuctuation edge of the spectrum. 

Below the spectrum of — ^A, the density of surface states is positive; hence, it is a Borel measure 
and the integrated density of surface states Ns{E) can be defined as its distribution function, i.e. Ns{E) = 
ns((— oo, E)) for E < —d. We will prove Lifshitz type behavior for Ns{E) for E \ Eq which is characteristic 
for fluctuation edges. However, the Lifshitz exponent, in the homogeneous case typically equal to — ^, is 
given by — ^ in our case. More precisely, we will show 

j.^ ln|ln(jV,(^))| _ di 
E\Eo ln{E -Eo) 2 ■ 
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1. The main results 



Let us now describe the general model we consider. Let H he a translational invariant Jacobi matrix 
with exponential off-diagonal decay that is H = ,^/g2d such that, 

(HO. a): h-j = hj for 7 G Z*^ and for some 7 / 0, /i^ / 0. 
(HO.b): There exists c > such that, for 7 G Z"', 

\hJ < -e^^lTl. 
c 

The infinite matrix H defines a bounded self-adjoint operator on Using the Fourier transform, it is 

easily seen that H is unitarily equivalent to the multiplication by the function 6 h{9) defined by 

h{e) = h^e'^^ where = {Oi, . . . ,60), 



acting as an operator on L^(T'^) where T*^ = /{2t:7/) (the Lebesgue measure on T'^ is normalized so that 
the constant function 1 has norm 1). The function h is real analytic on T"^. We normalize it so that it be 
non-negative and be its minimum. 

As both ends of the spectrum of our operator play symmetric parts, we only study what happens at a left 
edge, i.e. near the bottom of the spectrum. All our assumptions will reflect this fact. 

1.1. The case of a constant surface potentiaL We will start with a study of the density of surface 
states when the surfacic potential V is constant, i.e. V = fHg. We define the operator Ht = H + tl ® IIq. 
We prove two results on Ht- The first one is a criterion for the positivity of Ht and a description of its 
infimum when it is negative; the other result describes the density of the density of surface states near 
when Ht is non- negative. 
In the present section, we assume 

(HI): the function h : T'^ — > M admits a unique minimum; it is quadratic non-degenerate. 

If H is — ^A, then h = ho where 

(1.1) ho{d) := cos(0i) + • • • + cosiOd). 

In this case, assumption (HI) is satisfied. Below, we give an example why considering more general Hamil- 
tonians can be of interest. 

For the sake of definiteness, we assume the minimum of h to be 0. This amounts to adding a constant to H. 

We start with a characterization of the infimum of the spectrum of Ht- Therefore, write h{9) = h{9i, 62) 
where 6 = {61,62), 61 G T^S 62 G T^^^, define 

(1.2) I{6i,z)= [ T77r4^ d62. 



h{6i,62) - z 

We recall that the measures (162 is normalized so that the measure of ¥^^2 be equal to 1. 
We prove 

Proposition 1.1. Assume (HO) and (HI) are satisfied. 
Ht is non negative if and only if t satisfies 

(1.3) 1 +tIoo > where loo := sup / , . d62 

Assume now that 1 -|- t/oo < 0. Then, there exists a unique Eq G (-co, 0] such that 

y6ieT'^\ l + tI{ei,Eo)>0 and36ieT'^\ 1 + tl{6i, Eq) = 0. 

Moreover, Eq is the infimum of the spectrum of Ht . 

Proposition 1.1 is proved in section 5. 

Criterion (1.3) immediately gives the obvious fact that if t > then Ht is non-negative. As we assumed 
that h has only non degenerate minima, \i d2 = 1,2 and t < 0, then Ht is not non- negative. 
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We now turn to our second result. It describes the asymptotics of near when (1.3) is satisfied. 
Recall that iV* is the density of surface states of Hf. 



Theorem 1.1. Assume t satisfies condition (1.3). Define 

1 



h{0, 



One has 

• ifd2 = 1. 



•^0 di{di + 2)(2^)'^i JoetiQi - RQ^^R*) 



• ifd2 = 2: 

,.rt, . 2Vol{^'^^-^) £;i+<ii/2 
/ dNl{e) ~ ^ — — 

•^0 ^^^^ di{di + 2){2T:Y^^Det{Qi-RQ:2^R*) 

If d2>2) and 1 + t • / > 0, then, one has 

7o '^^£-0+ d{2^Y^/mtQ l + tl ^ ' 

If d2 > 3 and 1 + t ■ I = 0, if we assume, moreover, that 6i I(9i,0) := / {h{6i,92))~^d62 has a local 
maximum for 6i = 0, then one has 



if d2 = 3; 



E / Arg{-i\l-el\"^ + ~g{e^))dei 

di{di + 2)7r(27r)'^i J Det{Qi - RQ^^R*) 



• ifd2 = 4; 



^ r.rtr ^ 2V0/(S'^l-l) 

dNl{e) ~ : ^ — — 

di(di + 2)(27r)^i JZ)ei(Qi-iiQ^ii?*) 



i/d2 > 5; 



^i?^o+ d{2T:Y^/mtQ J ^ ' 



Here, we used the following notations: 

• Arg{-) denotes the principal determination of the argument of a complex number, 

• for n € {di, ^2}, is the n — 1 dimensional unit sphere, 

• g is a linear form defined below, 

• the function s and the constants c{di,d2) and J are defined by 

six) = l\x\'-^, c{di,d2)= f\'''~Hl-r^)^''^~^y'dr, J=[ titt!^^^! 

/ R*\ 

• Q is the Hessian matrix of h at that can be decomposed as Q = ( o ^ I • 

\R Q2J 



About the function g, it is defined as follows. We assume ^2 ^ 3 and 1 + tI = 0. Let /i2(^i) = inf h{9i,02). 

In section 5.1, we show that the function 9i ^ {h{6i, 62) — h2{0i))~^ d02 is real analytic in a neighborhood 
of 0. Using the Taylor expansion of this function near 0, one obtains 

^ + * / hfft h (ft . de2 = tg{ei) + om^). 

J 1^2 h{0i,d2) - h2{di) 
This defines the linear form g uniquely. Then, g is defined by 

m := (27r)'^VDet (Q2)5((Qi - RQ2^ R^T^'Hi). 

If the variables (^i, 62) separate in h, i.e., if 62) = /ii(^i) + /i2(^2)) the function g is identically 0. 

1.2. The case of a random surface potential. Let Vj^ be a random potential concentrated on the 
sub-lattice x {0} C [di is chosen as in section 0) of the form 

(1.4) K.(7i,72) = l"''' ^ n '7 = (71,72) G I^'^ x Z^^ = 

[0 if 72 7^ 0. 

and (<^7i)^^gz'*i ^ family of i.i.d. bounded, non constant random variables. 

Let uj± be respectively the maximum and minimum of the random variables (ti;^^)^^g^di , and let uj be its 
expectation. 

Finally, we define the random surfacic model by 

(1.5) H^ = H + K., 
and its IDSS by 

(n„ if) = E(tr(ni[v9(i?^) - v?(i/)]ni)) 
Following section 0, one regularizes into A'^, as in (0.5). 

Remark 1.1. An interesting case which can be brought back to a Hamiltonian of the form (1.5) with H 
and as above is the following. 

Consider F, a sub-lattice of 1/ obtained in the following way F = G({0} x 7/^) where G is a matrix in 
GSLf;(Z), the d-dimensional special linear group over Z, i.e. the multiplicative group of invertible matrices 
with coefficients in Z and unit determinant. One easily shows that the random operator 

if^(F) = -iA + J]a;^n^ 

(where 11^ is the projector onto the vector 5-^ S £^(Z^)) is unitarily equivalent to + where is defined 
in (1.4) and h{0) = hQ{G' ■ 0); here, /iq is defined in (1.1) and G' is the inverse of the transpose of G, i.e. 
G' = 

Definition 1.1. We say that E, an edge (or boundary) of the spectrum of Hlj, is stable if it is an edge of 
the spectrum oi H + tV^ for all t G [0, 1] . If an edge is not stable, we call it a fluctuation edge. 

Note that in the case of the introduction, this definition is equivalent to that given there. 
As in the introduction, one has to distinguish between 

(1) stable boundaries : at these boundaries, the IDSS is given by the IDSS of a model operator computed 
from the random model. 

(2) fiuctuation boundaries: at these boundaries, one has standard Lifshitz tails. 

To complete this section, let us give a very simple description of the spectrum of Hi^. One has 

Proposition 1.2. Let if^j be defined as above. Then 

a{H^) = IJ a{Ht). 

t£supp{Po) 

Here and in the following Pq denotes the common distribution of the random variables (a;-^2)72- 
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1.3. The stable boundaries. The stable boundary we are studying is the lower boundary that we assumed 
to be 0. Let us first give a criterion for the lower edge of the spectrum of H (that we assume to be equal to 
0) to be a stable edge. We prove 

Proposition 1.3. Write h{e) = h{ei,92) where 6 = {61,62), 61 G T'^S 62 G T''-. Then, is a stable spectral 
edge if and only if uj^ satisfies condition (1.3). 

Proposition 1.3 is an immediate consequence of Proposition 1.1 and Proposition 1.2. It gives the obvious 
fact that, if uj^ > 0, then is a stable edge. As we assumed that h has only non degenerate minima, we see 
that if ^2 = 1)2 and w„ < 0, then is never a stable edge. Actually, it need not be an edge of the spectrum 

Using the same notations as above, we prove 

Theorem 1.2. Assume (HO) and (HI) are verified. Assume, moreover, that is a stable spectral edge for 
H^. Then, one has 

(1.6) ifuj > 0, then liminf — J- > 1 and if ZJ < 0, then lim sup — i- < 1 

where is the IDSS of the operator with constant surface potential uJ, the common expectation value of 
the random variables (to-f-^ )'y^ . 

This result admits an immediate corollary 

Theorem 1.3. Assume (HO) and (HI) hold. Assume, moreover, that is a stable spectral edge for H^. 
Then, 



diidi + 2){2ttY^ J Det{Qi - RQ^^R*) 



ifd2 = l: 



• ifd2 = 2; 

Ns{E) ~ ^ ' — -. 

di{di+2){2nY^^Det{Qi-RQ:^^R*) 

Theorem 1.3 is an immediate consequence of Theorem 1.2 and the bound 

N'^-iE) < Ns{E) < N^+{E). 

As noted in the introduction, Theorem 1.2 is only necessary when = (in which case cJ > 0). Moreover, 
one obtains the analogue of (0.7) in the present case for ^2 ^ 3. 
The above results may lead to the belief that 

Ns{E) ~ NfiE) 

for all dimensions d2- Let us now explain why this result, if true, is not obtained for dimension d2 > 3. 
Therefore, we explain the heuristics behind the proof of Theorem 1.2; it is very similar to that of standard 
Lifshitz tails with one big difference when d2 > 3. 

Restrict H^ to some large cube. One wants to estimate the IDSS for H^^; for this restriction, this comes up 
to estimating the differences between the integrated density of states (the usual one) of the operator H^ and 
the integrated density of the operator Hi^_ (see Lemma 2.2). So we want to count the eigenvalues of H^ 
below energy E, say, subtract the number of eigenvalues of Huj_ below energy E, divide by the volume of 
the cube, and see how this behaves when E gets small. Assume ip is a normalized eigenfunction associated 
to an eigenvalue of ff^ below E. Then, one has {{H + Va;)(p,if) > E. Assume for a moment that is 
non negative. Then, we see that one must have both {Hip, ip) > E and ^p) > E. The first of these 

conditions guarantees that 99 is localized in momentum. So it has to be extended in space. If one plugs this 
information into the second condition, one sees that {Vi^ip, 93) ~ uJ with a large probability. So that, to (p, 
H^ roughly looks like H + uJIIq. There is one problem with this reasoning: 

as only lives on a hyper-surface, and as ip is flat, it only sees a very small part of ip; a simple calculation 
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shows that ||nQy3|| ~ ^i^^ other hand, when one says that ip is roughly constant, one makes an 

error of size E°' (for some < a < 1); hence, for dimension d2 > 3, this error is much larger than the term 
we want to estimate, namely, {V;^ip,(p). In other words, because (p is very flat, we can modify it on the 
hyper-surface (e.g. localize the part of it living on the hyper-surface) with almost no change to the total 
energy of (p; hence, we cannot guarantee that is also flat on the hyper-surface, which implies that {Vuj'P, f) 
need not be close tJ with a large probability. 

1.4. The fluctuation boundaries. In this section we assume that the infimum of E which we call Eq is 
(strictly) below mf{a{H)), so that Eq is a fluctuation edge. In this case, we consider a "reduced" operator 
H which acts on i'^iZ'^'^). In Fourier representation this operator is multiplication by the function h given 
by: 

(1-7) HOi) = ([ T7^4r—Erde2) + E, 



\JTd2 h(9i,62) - Eq J 

We will reduce the proof of Lifshitz tails for = if -|- to a proof of Lifshitz tails for the reduced operator 
= H + Vuj (where is a diagonal matrix with entries {uj-^j^)-^^). To prove Lifshitz tail behavior for H^j 
we have to impose a condition on the behavior of h near its minimum. We either suppose: 

(H2): the function h : T'^i ^ R admits a unique quadratic minimum. 

or we assume the weaker hypothesis: 

(H2'): the function h : T"^ — > M is not constant. 

Moreover, we always assume that the random variables defining the potential (0.1) are independent with 
a common distribution Pq. We set w_ = inf(supp(Po)) ^ind assume: 

(H3): Pq is not concentrated in a single point and Pq{[oj^,uj^ + e)) > C for some k. 

We will prove below: 

Theorem 1.4. // (H2) and (H3) are satisfied then 

ln|ln(iV,(^))| _ di 
E\Eo ln{E - Eo) 2 ■ 

We have an additional result for low dimension of the surface: 
Theorem 1.5. Assume (H2') and (H3) hold. If di = 1 then 

^. ln|ln(iV,(^))| ln{n{E)) 
hm — — 7— — -- — = — lim — -— 

E\Eo ln{E - Eo) E\Eo {E - Eq) 

where n[E) is the integrated density of states for H. 
If d2 = 2, then 

^. ln|ln(iV,(S))| 

lim ^ — = —a 

E\Eo ln{E - Eq) 

where the computation of a is explained below. 

For the sake of simplicity, let us assume Eq = 0. The Lifshitz exponent a will depend on the way h vanishes 
at S = {9i\h = 0} and on the curvature of S. 

To describe it precisely, we need to introduce some objects from analytic geometry (see [19] for more 
details). If i5 is a set contained in the closed first quadrant in then its exterior convex hull is the convex 
hull of the union of the rectangles Rxy = [x,oo) x [y,oo), where the union is taken over all {x,y) G £. 

Pick 9q & S and consider the Newton diagram of h at Oq, i.e., 

(1) Express h as a Taylor series at ^o, HO^, 0^) = Eij aij{9^ - 9^y{e^ - 9^)^, 9 = {6^,9^). 

(2) Form the exterior convex hull of the points {i,j) with Oij 7^ 0. This is a convex polygon, called the 
Newton polygon. 

(3) The boundary of the polygon is the Newton diagram. 



The Newton decay exponent is then defined as follows. The Newton diagram consists of certain line segments. 
Extend each to a complete line and intersect it with the diagonal line 9^ = 6^. This gives a collection of points 
(ofc, Ofc), one for each boundary segment. Take the reciprocal of the largest and call this number a{h, 6q); 
it is the Newton decay exponent. Define a{h,9o) = min{d!(/i o To, ^o) : ^o(') = 9o + T{- — 9o), T G 5*^(2, M)}. 

Similarly, define a{h, 0) if 6 is any other point in 5, the zero set of h. Then, the Lifshitz exponent a is 
defined by 

(1.8) a = mma{h,0). 

The Lifshitz exponent a is positive as i— > a{h, 0) is a positive, lower semi-continuous function and S is 
compact (see [19]). 

Remark 1.2. Let us return to the example given in Remark 1.1. In the section 6, we check that (H.2') 
holds in this case; so for d = di + d2 = 2>, Theorem 1.5 applies. 



2. Approximating the IDSS 

To approximate the IDSS, we use a method that has proved useful to approximate the density of 
states of random Schrodinger operators, the periodic approximations. We shall show that the IDSS is well 
approximated by the suitably normalized density of states of a well chosen periodic operator. 

2.1. Periodic approximations. Let (w7i)^jggdi be a realization of the random variables defined above. 
Fix G N*. We define ^ a periodic operator acting on £"^{7/) by 

/32e(2iv+i) ■ 



Here, U^^^j^^ = 11^ /{2N + 1)^*^, 5i = {^ji)j^^d is a vector in the canonical basis of ^^(Z'^) where 6ji is the 

Kronecker symbol and, d = di or d = d2, the choice being clear from the context. As usual, is the 

orthogonal projection on a unit vector u. 

By definition, is periodic with respect to the (non degener ate) lattice (2^ + 1)^^^. We define the density 
of states denoted by as usual for periodic operators: for ip G C^(M), 

((/p,dn^) = / (p{x)dn^{x) = lim . ^ V {6^,(p{H^)6^). 

h\<L 

This limit exists (see e.g. [4, 23]). In a similar way, one can define the density of states of H; we denote it by 
dnQ. The operators {H^)uj,n are uniformly bounded; hence, their spectra are contained in a fixed compact 
set, say C. This set also contains the spectrum of H^^ and H. We prove 

Lemma 2.1. Pick U CE. a relatively compact open set such that C CU. There exists C > 1 such that, for 
ip £ C^(IR), for K gN, K > 1, and N £ N*, we have 



(2.1) {ip, dn) - {2N + 1)'^^E{((^, [dn^ - dno])} 



[CK\ 



K 



0<J<K+d+2 



dJ 



Proof of Lemma 2.1 Fix ip G C^(M). As the spectra of the operators are contained in U, we may 
restrict ourselves to (p supported in U which we do from now on. By the definition (0.2), one has 

(2.2) {^,ns) = E ^ (5o ® 572> IfiHu.) " v{H)]6o S^,) = Mn{^) + Rn{v) 



where 



\|72|<Af 



\l72|>iV 



Let us now show that 
(2.3) 



\Rn{v^)\ < 



CK \ 



K 



sup 

0<J<K+d+2 



Therefore, we use some ideas from the proof of Lemma 1.1 in [17]. Helffer-Sjostrand's formula ([10]) reads 

ip{H^) = ^ f ^{z) ■ {z - H^r^dz A dz. 

^71" Jc 

where (f is an almost analytic extension of (p (see [22]), i.e. a function satisfying 

(1) for z G M, (f{z) = ip{z); 

(2) supp((^) C {z G C; |Im(z)| < 1}; 

(3) ^gSHzGC; |Im(z)| <1}); 

d ~ 

(4) the family of functions x 1— > -g^{x + iy) ■ \y\~'^ (for < \y\ < 1) is bounded in 5(]R) for any n G N; 
more precisely, there exists C > 1 such that, for all p,q,r G N, there exists Cp^q > such that 



(2.4) 



sup sup 

0<\y\<l xm 



dz 



{x + iy) 



< C'Cp^q sup sup 
q'<r+q+2xeM 
p'<p 



As we are working with (p with compact support in U, its almost analytic extension can be taken to have 
support in {U + [—1, 1]) + z[— 1, 1] (see e.g. [()]). 

We estimate E(|((5o S-y^, [(p{Huj) — (p{H)]6o ^72)!) for I72I > N. Using the fact that the random variables 
('^72)72 bounded, we get 

K{\{5o 5^,, [ip{H^) - ip{H)]5o ® <572)l) 



< — E 

71 eZ'*! 



dz 
dip 
dz 



|(5o ® ^72, {{z - H^)-^ -{z- H)-^) So 6-y,)\dxdy 
■E{\{do(g)d^^,{z - H^y'^d^^ (g)6o)\ ■ |((5^i ®6o,{z- H)-^5o®6^^)\)dxdy 



where z = x + iy. 

By a Combes-Thomas argument (see e.g. [18]), we know that there exists C > 1 such that, uniformly in 
(u!^)^, 71 G Z*^! and > 1, we have, for Im(z) 7^ 0, 



(2.5) \{6^, 5^„iz - H^r^5^> 'S)6yJ\ + \{5^, ® 6^„iz - H)-'6y^ < 



"7i ^ "72 ' 



< 



c 

\Im{z) 



|Im{^)|{|7i-7;| + |72-7^|)/C 
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Hence, for some C > 1, 



\Rn{v)\<C 



dz 



< C 



dz 



\lm{z 



^^e-|Im(.)(|7i| + l72|)|/C^^^ 
|lm(2:)|^ 



Taking into account the properties of almost analytic extensions (2.4), for some C > 1, for K > 1 and 
> 1, we get 



\Rn{v)\ < C 



K+l 



< 



CK \ 



(w+[-i,i])+ihi,i] 



sup 

0<J<K+d+2 



y\Ke-\y^\/Cdxdy sup 
(x) 



0<J<K+d+2 



d-^if 



dJ 



X 



(x) 



X 



This completes the proof of (2.3). 

We now compare Mn{^) to {2N + iy^E{{ip, [dn^ - dn^])}. Therefore, we rewrite this last term as 
follows. Using the (2A^ + 1)1/ periodicity of and H, we get 



7ez<^ 

|7|<Ar+L(2Ar+l) 



7eZ'* 
|7|<Ar 



This gives 



(2.6) 



{2N + ir{y,,dn^)=E 



\|7|<7V 



On the other hand, as the random variables (^^2)^2 are i.i.d. and as H is Z'^-periodic, as in [18], one 
computes 



E 



/ 

\|7|<iV 



E 



\< 
< 



72/ 



7iGZ'*l, |7i|<Ar 
,72eZ'*2, |72|<Ar 



(2iV + l)'^iE 



72eZ'*2 

\l72|<Af 



Combining this with (2.6), we get 



/ 



{2N + lf^E[{^,dn!;^)]=E 



^72) 



, 72eZ''2 

\|72|<iV 
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Of course, such a formula also holds when is replaced with H. In view of (0.2), (2.3) and (2.2), to 
complete the proof of Lemma 2.1, we need only to prove 



(2.7) 



E 



72? 

h2\<N 



< 



CK 



K 



sup 

0<J<K+d+2 



(x) 



for if, K J and N as in Lemma 2.1. 

Proceeding as above, for 72 G Z'^^^ |^2| < A^, we estimate 



< C 



71 « 



.7^e({2Ar+l)Z'*2)* 7^=0 



dxdy- 



¥.[\{5o®5^,,{z-H'::r'6^,®5,,^)\- 



Here we used the fact that the operators and coincide in the cube {I7I < N}. 

As H^^ satisfies the same Combes-Thomas estimate (2.5) as , doing the same computations as in the 
estimate for R]\f{{p), we obtain (2.7). This completes the proof of Lemma 2.1. □ 

Obviously, one has an analogue of (2.1) for 

^s,norm; norm" One needs to replace and H with 

their obvious counterparts i.e., choose the random variables (^^2)72 to be the appropriate constant. 
This enables us to prove 

Lemma 2.2. Fix I, a compact interval. Pick a > 0. There exists vq> Q and p > such that, for 7 G [0, 1], 
E £ I , v £ {0, z^o) and N > u'P , one has 



< Ns{E) <E{N^ {E + u)) + e-" 



(2.8) 

where N^^^^^ = — N^_, and (resp. N^_) is the integrated density of states of (resp H^_, i.e. 
where ojy = uj- for all ^). 



N 



or that of 



Let us note here that one can prove a similar result for the approximation of A'^s^norm by 
Nl by iVi'^. 

Proof Let us now prove Lemma 2.2. Pick if a Gevrey class function of Gevrey exponent a > 1 (see [11]); 
assume, moreover, that ip has support in (—1, 1), that < < 1 and that = 1 on (—1/2, 1/2]. Let E £ I 
and 1/ G (0, l),and set 

Then, by Lemma 2.1 and the Gevrey estimates on the derivatives of (p, there exist C > 1 such that, for 
N > 1, k > 1 and < u < 1, we have 



* (j9 I - 



(2.9) 



We optimize the right hand side of (2.9) in k and get that, there exist C > 1 such that, for > 1 and 
< < 1, we have 



\E{{^E,u,dN!^)) - {^E,u,dNs)\ <C{N + v 're 



Now, there exist i^o > such that, for < z/ < z/q and N > u ^, we have 

tNw I - -TAT M ^ __!y-'7/(2a) 



l'|3p-(Afi'/C)i/(i+°)+C(Ari^/C)-i/(i+''' 
1-r, 



(2.10) 



\E{{^E,vAN':))-{^E,vANs)\<e 
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By definition, (fE.v = 1 on [0,ii^], and ^e,v has support in [—v,E + u] and is bounded by 1. As dN^ and 
dNs are positive measures, we have 

(2.11) m^A^)) < mvE,u, dN^)) < m!!JE + u)). 

Hence, by (2.10) and (2.11), we obtain 

Ns{E) < {ipE,,,dN,)=E[{^E,.,dN^)] + [{^E,.,dNs)-E{{ipE,u,dN^))] 
<E(iVi^(i? + ^)) + e-^''''''"^ 

and 

Ns{E) > {ipE-,,u,dNs) = E[{ipE-u,.,dN^)] + [{ipE-u,.,dN,) -K{{ipE-,,u,dNff))] 
>E(iV5(i?-i.))-e--~'"*'"^ 
Tliis completes tlie proof of Lemma 2.2. □ 

2.2. Some Floquet theory. To analyze the spectrum of , we use some Floquet theory that we develop 
now. We denote by J" : L'^{[ — 7r,7r]'^) — > (^{TL'^') the standard Fourier series transform. Then, we have, for 
t.GL2([-7r,7rf), 

{H^u){B) = {T*H^Tu){e) = h{e)u{e) + ujii{iiiiu){e) 

where iUpu){e) = j^e'f"' [ e~'^\{e)de. 
Define the unitary equivalence 

where the iU'y{6))y^^d^ are defined by 

(2.12) u{6) = e^'^^u.^{6) where the functions {6 ^ u-^/{6))^^^d_^^^ are ^^^^ — ^Z'^-periodic. 



The functions iu'y)^^^^^ ^ computed easily; if the Fourier coefficients of u are denoted by {u^)^(=,id, then, 
one gets 

(2.13) u,{e) = J2 V(2^+i)/3e^^'^+^)'^^ 

l3£Zd 

The operator UT*H^TU* acts on L'^{[-^2W+i^ IWTi^'^) ® ^^(^ 2Ar+i); it is the multiplication by the matrix 

(2.14) M!^{e)=H^{9) + V^ 
where 



(2.15) H (0) = ((/i^_^/(0)))(^^^,)g(2d^^^)2 and = ((^/3i^/3i/31'^/320'5/3^o))(^^^^/)g(zrf]^^^)2. 

Here, the functions {h^)^^id^ ^ are the components of h decomposed according to (2.12). The (2A^ + l)'^ x 

{2N + l)'^-matrices H^{9) and are non-negative matrices. 

This immediately tells us that the Floquet eigenvalues and eigenvectors of with Floquet quasi- 
momentum 9 (i.e. the vectors, u = {up)p^id), solution to the problem 

[h^u =Xu, 

\up+^ = e-^^'^up for /3 G Z'^, 7 G (27V + 1)1'^) 
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are the eigenvalues and eigenvectors (once extended quasi-periodically) of the (2A^-|- l)'^ x (2A^-|- 1)"^ matrix 
M^{e). For E eR, one has 

M^{E)= f dn^{E)= [ ^{eigenvalues of MJ(0) in [O,S]}d0. 

•'I 2N+l'2N+li 

Considering H as {2N + l)Z'^-per iodic on Z'^, we see that the Floquet eigenvalues of H (for the quasi- 
momentum 6) are {h{6 + 2^+1 ))^gz^^ ^5 Floquet eigenvalue h{6 + 2^+1) ^® associated to the Floquet 
eigenvector Uj{9), 7 G ^2Af+i defined by 

uJ9) = ^_(e-^(''+2?^)/3) . 

(2A^ + ^pG^2]v+i 

In the sequel, the vectors in £^(Z2jy_,„;^) are given by their components in the orthonormal basis {uy{9))^^.^d 
The vectors of the canonical basis denoted by {vi{0))i^'^d have the following components in this basis 

We define the vectors {vi)i^2,^^ by 

3. The proof of Theorem 1.2 
To prove Theorem 1.2, we will use Lemma 2.2 and the Floquet theory developed in 2.2. We will start 



with 



3.1. The Floquet theory for constant surface potential. We consider the operator = where 

71t^2JV+l 

form (2.14) where 



LO = (i) di is the constant vector and t ^ 0. The matrix {9) defined by (2.14) for takes the 



'2JV+1 

{i32A)^c4l+i? 



Our goal is to describe the eigenvalues and eigenfunctions of M^{9). As usual, we write 9 = (0i,02)- By 
definition, the operator ijf is Z'^^ x {2N + l)Z'^2.pei.iodic. It can be seen as acting on f{Z'^\e^{Z'^'')); as 
such, we can perform a Floquet analysis in the 0i-variable as in section 2.2 (in this case, just a discrete 
Fourier transform in ^i) to obtain that H^^ is unitarily equivalent to the direct sum over 61 in T^i of the 
2N + 1-periodic operator {9i) acting on ^^(Z'^^) defined by the matrix 

72G(2Af+l)Z<*2 

Here h(9i; (32) is the partial Fourier transform of h{9i,92) in the 02-variable. 

For each 9i, we now perform a Floquet reduction for {9i) to obtain that H^^ [Oi] is unitarily equivalent 
to the multiplication by the matrix 

M^{e,^92) = ((M^i,^2;/32 -/3^) + M,.o<5,^o))(,,,,,),(,^^^^)2 

The matrix-valued function {61,62) ^ M/^(^i,02) is 27rZ'^i -periodic in 9i and gj^q^x^'^^-periodic in ^2- 
It is a rank one perturbation of the matrix Mq {9i,92)', the eigenvalues of this matrix are the values 

h [9i, 92 + ) . Let us for a while order these values increasingly and call them {E^ (^1, ^2, i))i<n<n.jv) 



2N + l_ 

where tcat < (2A^ -|- 1)*^^ ^-^^g (\^q ^loi repeat the eigenvalues according to multiplicity). Assume t > 0. The 
standard theory of rank one perturbations [24] yields 

Lemma 3.1. Forl<n< un , if E^ {9i,92,0) is an eigenvalue of multiplicity k, then 
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• either it is an eigenvalue of multiplicity k for M^^ {9i); 

• or it is an eigenvalue of multiplicity k — 1 for (Oi) and the interval {E^ {9i,0), E^_^i{6i,92,0)) 
contains exactly one simple eigenvalue; this eigenvalue is given by the condition 

t{5o,iE-M^i6i)y^6o) = l; 

Here, we took the convention 02> 0) = +00. One has a symmetric statement for t < 0. 

For 1 < n < UN, let {^n ji^i^^'ij^))i<j<jn denote orthonormalized eigenvectors associated to the eigenvalue 
E^ {6i,02,t) where j„ denotes its multiplicity. 

In the sequel, it will be convenient to reindex the eigenvalues and eigenfunctions of the matrix {9i,92) 
as {E^(ei,e2,t)) ^„d2 and (^i, ^2, 0) ^1^2 ■ Clearly, the functions (^1,^2) ^ E^jei,e2,t) and 

72fc*^2JV+l 'y^^ 2N+1 

{61^62) I— > ipj^{9i,92,t) can be chosen to be 27rZ'^i -periodic in 9i and ^j^^i.^ ^^^ -periodic in 92- 
Let us now show the 

Lemma 3.2. The eigenvalues of M^^ {9) are the values {E^-^^-y^{9i,92,t); 71 G Z2]y_,_]^, 72 G where 
(3-2) £'^1,^2 (6*1, 6I2, t) = (^9i + ^^ '^2,^ 

A normalized eigenfunction associated to the eigenvalue i?^ + 2N+1 > ^2 ; ^ j is the vector 
(3.3) v,,,M,02,t) := (2iV + 1)-^^/^ (e-^^^'^^^^^^g (9, + ^,^2,t)) ^ , 

i.e. the vector of components 

(3.4, 1)-*/= (,-«.<»..|a.^. + 2|tt))«4V„ 

if (9\,t) has components {(^l(9\)) • 

T/ie vectors (^^1,72(^1) ^2, /orm an orthonormal basis of £'^{7j2n_^_^ x TJ^j^^^-^). 

72eZ2Ar+l 

Proof. Orthonormality is easily checked using the fact that the vectors {(p^ {9i,92,t)) 0.2 form an 
orthonormal basis. 

Let us now check that u^^^^jC^i; ^2) satisfies the eigenvalue equation for M^^ [9) and E^j^^y^{9i,92,t) given 
in (3.2). Therefore, first note that the matrix M^^ (9) is nothing but the multiplication operator by the 
matrix- valued function {9i) to which one has applied the Floquet reduction of in the ^i-variable. Hence, 
by (2.13), the matrix elements of {9) given by (2.15) satisfy, for j3i G Zj^y+i, 



(3.5) 



Both sides in this equality are matrices acting on ^^(Zg^^-^), the matrices rnp^_pi^{9) being defined as 

mp,-p'^{9) = ((V-/3l,/32-/3^(^)))(^2,/3^)G(<^+i)2- 
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If we now apply both sides of equation (3.5) to the vector (^i + 2N+1 ' ^ ) ' obtain, for /3i E ^2Af+i 



* V 2iv + 1 y V 2A^ + 1 ' 

This is nothing but to write 

This completes the proof of Lemma 3.2. □ 
In the course of the proof of Theorem 1.2, we will use the 

Lemma 3.3. Fix t such that t > if d2 = 1,2 and 1 + t/oo > i/ ^2 > 3. Then, for p > 2, there exists 
C > such that, for N > E^p and E sufficiently small, the eigenvalues of M^^ satisfy 

(3.6) E^,,,,{ei,e2,t)<E^(]-^\ <CE 



, 2iV + 1 , 

Proof. When t is positive, (3.6) is clear by Lemmas 3.2 and 3.1, that is, by the intertwining of the eigenvalues 

of M^{e) and M/^((9), and as the eigenvalues of M^{e) are the values h (^Oi + ^^,^2 + 2^+1 ) ^^^^^ 

satisfy (3.6) as h{e) > C\e\'^. 

Assume now that d2 > 3 and t satisfies 1 + t/oo > 0. To complete the proof of Lemma 3.3, by 
Lemma 3.2, it is then enough to prove that, there exists C > such that for, one has 

|ei|2>c^^V72, E^^{6i,e2,t) > E. 

By the intertwining properties and the properties of h, this is clear except for the lowest of the {E^^)^^. 
Assume now that |(9i|2 > E. Then, by our assumptions on the behavior of h near its minimum, for some 
C > 0, one has that (^1, e) 1-^ I{Oi, e) is real analytic in > E} x {|e| < E/C}. Hence, using a standard 

estimate for Riemann sums, we get that, for |6'i|2 > E and |e| < E/C, 

1 + t{5o, {M^iOi) - ey^do) = 1 + tl{ei,e) + 0{E~^EP) 

So that, as 1 + t/oo > 0, for E sufficiently small, the equation l + t{5o, {M^ {9i) — e)~^ 60) = has no solution 
for |6'ip > E and |e| < E/C. By the above discussion, this implies that, all the E!^^{9\,92,t) lie above E/C. 
This completes the proof of Lemma 3.3. □ 

3.2. The proof of Theorem 1.2. We now have all the tools necessary to prove Theorem 1.2. Notice that, 
as tJ > u^, as 1 + oj-Ioo > 0, we know that 1 + ujI^ > 0. So that the asymptotics for Nf^E) are given by 

Nf{E)^r^^^C{uJ)-f{E). 

The precise value of the constant C{io) and of the function f{E) are given in Theorem 1.1. The constant 
C(uJ) is a continuous function of tJ; and, for any c G M, the function f{E) satisfies f{E + cE'^) ~ f{E) when 
E ^ moreover, / is at most polynomially small in E. All these facts will be useful. 

We start with the proof of (1.6). We will use Lemma 2.2. As above, fix N large but not too large, say 
_/V ~ E~P for some large p. Fix 5 > small. Consider the matrix M^j^^{9) obtained by the Floquet reduction 
of + (a7 + (5)nQ. Let Ti^ {E,9) be the spectral space of M^^(6') associated the eigenvalues less that E. 
Then, we prove that 
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Lemma 3.4. Fix 5 > 0, p > 2 and a G (0, 1/2). For N ^ E ^ and E sufficiently small, with a probability 
at least 1 — e^ for all 6 and all ip G 7ig^ {E, 0), one has 

{M!^{e)^,ip)<E\M\ 

This lemma immediately imphes the desired lower bound. Indeed, it implies that, for N ~ E~p, with a 
probability at least 1 — e^ " , one has 

N^+5iE) = [ tt{eigenvalues of M^+^C^) in [0, E]}d6 

h TT TV Td 

I 2JV + 1 ' 2JV + 1 J 

< / tj{eigenvalues ofM^{9) in [0,E]}de 

I 2JV+1 '2JV+lJ 

Taking the expectation of both side, and using (2.8) (and the fact that the number of eigenvalues of M^{6) 
and M^^g{e) are bounded by (2A^ + 1)^), we obtain 

Nf+\E -E^)- CE'^Pe-^"' < N^iE) 
Considering the remarks made above, we obtain 

C(uJ) <liminf^^^. 
^ ^ - E^o+ f{E) 

As C{uJ) has the same sign as ZJ, this completes the proof of (1.6). 

Proof of Lemma 3.4. Pick E small and ip G Tif {E, 9). Then, by Lemma 3.3, ip> can be expanded as 

p = 'y ^ '^71,72^71,72 (^) ^) 

|7i|<C£;i/2Ar 

rT/dn 

726^2^^^ 

where the vectors {v^{6))^ are given by (3.3) and (3.4). Using these equations, we compute 
(3.7) (K.V¥')= E '^ftl^/^il' 

where 

(^•^^ -^/3i = f2/V+\yi/2 E e^^c^i and c^, = ^ «7i,72 60, < (^1 + ^2, t) 

So the vector (A^J^^ is the discrete Fourier transform of the vector c = (c^i)-yi supported in a ball of radius 
CE^/'^N . To estimate this Fourier transform, we used the following result 

Lemma 3.5 ([18]). Assume N , L, K, K' L' are positive integers such that 

• 2iV + 1 = {2K + 1)(2L + 1) = {2K' + 1)(2L' + 1) 

• K < K' and L' < L. 

Pick a = (a„)„g^d G ^^(Zg^ , -^) such that, 



for \n\ > K, a„ = 0. 



Then, there exists a G i'^{'Ij2^_^j) such that 



(1) \\a - a||^2(zd^^^) < CK,K'\\a\\i2(^zl^^^) where Ck,k' ^k/k'~.o K/K'. 

(2) write a = {aj)j^xd^ ^; for I' G Z2L'+i '^'^^ ^' ^ ^2A"+i' have 



X - , 2^j.(i' + fc'(2£' + l)) , ^ 2nrk' 

y ^ djC 2iV+l = 2^ CLjC^K^, 
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(3) M\p{ii^^^) - Mpm^^.y 

This lemma is a quantitative version of the Uncertainty Principle; it says that, if a vector is localized in a 
small neighborhood of (here, of size K/N), up to a small error 6, its Fourier transform is constant over 
cube of size N/{6K). 

To apply Lemma 3.5, we pick N such that {2N + 1) = {2K + 1) (2L' + 1) (2M + 1) where K > CE^/'^N; 
this is possible as ~ E~p with p large; we pick for example, L' ~ CE~^^~'^^/'^ and L' ~ CE~'^/'^ (for some 
fixed < 1/ < 1). 

We apply Lemma 3.5 to the vector c = (0^1)71 defined in (3.8); by Lemma 3.5, there exists c = (£7^)71 so 
that, if we set 



e 2iv+i c. 



■71 



TiG^ajv+i 



then, for 7I G ^2L'+i ^'^'^ /'i ^ ^2/^'+i' have 
(3.9) 



^7(+/3i(2L'+l) - -4^i(2L'+l)- 

Fix r/ > small to be chosen later. We replace A by ^4 in (3.7) and use the boundedness of the random 
variables to obtain 

Using (3.9) and points (1) and (3) of Lemma 3.5, we get that 

-^^/2 + \ ^ 

7] {2U + lY^ 



c 



\ 



{2L' + ir^\Af,,{2L' + l)y 



Pick rj such that t] ■ u;+ < 5/4: and E sufficiently small that ^E^/"^ < 5/4. We then obtain 



(3.10) (K.V^)< E 
Now, if uj satisfies 



6/2 + 



1 



'2L' + l)'^i 



{2L' + lf^Ap,{2L' + l) 



(2L' + 



J y 

/ ^ l)di 



< 07 + 5/2 



then, (3.10) gives 



(y>,y.)<(^ + 5) E (2L' + l)'^H^/3i(2i' + l)l' = (^^+5y',^) 



2if' + l 



where is defined in (3.1). Here, we have used the points (2) and (3) of Lemma 3.5, and the definition (3.8) 

of the vector c = (07^7^. 

Summing all this up, we have proved 

Lemma 3.6. Pick < < 1. Pick N as described above. For E sufficiently small, the probability that, for 
all 9 and all if G Ti^{E,6), one has 



{Mf^{e)ip,ip)<E\\ip\\\ 



is larger than the probability of the set 



ydl 



^2K'+V (2L' + l)'ii 



E ^l[+f3[i2L'+l) < ^ + 6/2 
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The probability of this event is estimated by the usual large deviation estimates (see e.g. [7, 5]). This 
completes the proof of Lemma 3.4. □ 



4. The fluctuating edges 

In this section, we investigate the behavior of the density of surface states Ns{E) at the bottom Eq of 
the spectrum of H^^ in the case when Eq < inf (t{H) = 0. As we saw in Section 0.1, this is always the case 
for dimension d2 = 1 or d2 = 2 and it holds in arbitrary dimensions if the support of common distribution 
Pq of the has a sufficiently negative part. Thus, we are looking at a fluctuation edge as described in 
Section 0.2. Due to the symmetry of the problem we may, of course, consider the top of the spectrum in an 
analogous way. 

4.1. A reduced Hamiltonian. In the present situation it is convenient to think of the Hilbert space 
^2^^di+rf2) as a direct some of ^^(Z'^i x {0}) =: Tib and ^2(zrfi+<i2 \ ^di x {0}) =: Tis, the indices referring 
to "bulk" and "surface" respectively (see [13] whose notations we follow). According to the decomposition 
TC = Tis ffi T~(-b we can write any operator A on H as a matrix 

[ Ahb 

where Agg and A^b act on Tig and Tib respectively and Agb : Ti-b Tis, Abg : Tig Tib "connect" the two 
Hilbert spaces Tig and Tib- The bounded operator A is symmetric if A^g = Agg, A^^^ = Abb and A*^ = Abg. In 
the case of our random Hamiltonian we have: (-ffa;)ss = (-ffo)ss + while {H^)bb = {Ho)bb and Hgb as 
well as Hbg are independent of the randomness. Moreover, by assumption, {Hi^)bb > 0, while mia{H^) < 0. 
Consequently, the operator ((-ffo)66 — ElLbb)~^ exists for all £' < and the operator 

Gg{E) := {Ho)gg + K. — Hgb {{Ho)bb — Elbb) ^ Hbg — Eigg 

the so called resonance function is well defined. The operator Gg{E) is a sort of a reduced Hamiltonian. Its 
inverse plays the role of a resolvent. It is not hard to show that the set 7?.(ff(^) = {E s] — oo,0[;0 E a{Gg{E))} 
(the resonant spectrum) agrees with the negative part of a{H^). See Prop. 1.2 in [13] for details. For later 
reference, we state this as a lemma: 

Lemma 4.1. For E < 0, E is an eigenvalue of if and only if is an eigenvalue at Gg{E). Moreover 
the multiplicities agree. 

In fact a little linear algebra proves that, for block matrices, we have 

, . ( A B\~^_( (A-BD-^C)-^ -A-^B{D-CA'^By^ \ 

^ ^ \C DJ ~ \ -D-^C{A + BD-^C)~^ (D-CA-^By^ J 

when all the terms make sense. 

We denote by N{A, E) the number of eigenvalues (counted according to multiplicity) of the operator A 
below E. For Al = [—L,LY we set {H^^L)ij = {H^)ij if i,j G A^ and {H^^L)ij = otherwise. For energies 
E below zero the integrated density of surface states of is given by 

Defining 

Gs{E) = {H^x)ss — {HL)sb{{HL)bb - Elbb) ^{HL)bs — E'^gg. 

We have, as above, that < is an eigenvalue of H^^ l if and only if is an eigenvalue of Gg{E). 

In the following, we will express the density of surface states Ng{E) (for £^ < 0) in terms of the operators 

gHe). 

Lemma 4.2. The eigenvalues Pn{E) of Gg{E) are continuous and decreasing functions of E (for E < 0). 
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Proof of Lemma 4.2. Continuity is obvious from the explicit form of the entries of the (finite-dimensional) 
matrix Gg{E). In the following calculations we omit the superscript L. Let > E2 > Ei, then 

Gs{Ei) — Gs{E2) = —Hsbi{Hbb — El) ^ — {Hbb — E2) ^)Hbs — {Ei — E2) 

= {E2 — Ei)Hsb{{Hbb — El) ^{Hbb — E2) ^)Hbs + {E2 — El). 

Since Ei,E2 < the operator {Hbb — Ei)~^{Hbb — -^2)"^ is positive, the operator Gs{Ei) — Gs{E2) is 
positive. □ 

Proposition 4.1. For E < 0: 

N{H^^L,E) = N{G^{E),0). 

Proof of Proposition 4.1. For E sufficiently negative, Gs{E) is a positive operator. Let us now increase 
E (toward £" = 0). Then, E is an eigenvalue of H^^^l if one of the eigenvalues of G^{E) passes through zero 
and becomes negative. □ 



Ns{E) = lim iV(Gf(^),0) 

L— >oo 



It follows from this proposition that (for E <{)) 

Gg{E) depends on £^ in a rather complicated way through the resonance function. We will therefore 
approximate G^{E) by an operator with much simpler dependence on E in the following way: 
let Eq = inf a{H^) then we set: 

GsiE) = iH^^L)ss - {HL)sb{{HL)bb - Eo)~^{HL)bs - E 
This operator should give a good estimate for the eigenvalues of H^^ near Eq, in fact: 

Lemma 4.3. For Eq < E < : 

N{Gs{E),0) < N{Gs{E),0) 

Proof of Lemma 4.3. 

Gs{E) — Gs{E) = {E — Eo)Hsb {{Hbb — E) ^{Hbb — Eq) ^) Hbs- 

So 

G^E) > G^{E). 



□ 



For a bound in the other direction we observe that: 

Lemma 4.4. For Eq < E < Ei < we have 

Gs{E)-Gs{E)<C{E-Eo) 

Remark: The constant C in the above estimate depends on Eq and Ei. 
Proof of Lemma 4.4. 

Gs{E) — Gs{E) = {E — EQ)Hsb{{Hbb — E)^^{Hbb — EQ)~^)Hbs 

< {E — EQ)Hsb{{Hbb — Ei)~^{Hbb - EQ)~^)Hbs 
<G{E-Eq). 

Here, we used that 

{Hbb-E)-^ < {Hbb-Ei)~\ 

Summarizing, we have got: 

Proposition 4.2. There is a constant C , such that for Eq < E < Eq/2 < 

N{Gs{E),0) < N{H^,E) < N{Gs{E) - G{E - Eq),0). 
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The advantage of having Gs{E) rather than Gs{E) lies in the fact that Gs{E) depends hnearly on E, in 
fact: 

Gs{E) = Hss — HsbiHbb — Eq)~^Hi)s + Vuj — E 
= H + V^-E 

where H is the operator 

H = Hss — HsbiHhh — Eo) ^Hbs- 
This operator is of a similar form as the Hamiltonian H, however it acts on ^^(Z"^!), i.e. on the surface only 
where the random potential lives. The price to pay is the complicated looking "bulk term" HsbiH^b — 
Eo)~ Hbs- 

Nevertheless, H is still a Toeplitz operator and it is not too hard to compute its symbol, i.e. its Fourier 
representation. 

In fact, a look at formula (4.1) shows that 

(4.2) H=[{[H-E^)~^)J-' + Eo. 

Consequently the symbol of H is given by: 

We summarize these results in a theorem: 

Theorem 4.1. Let H^; = H + Vuj as in (1.5) satisfying assumption (HI). Assume moreover, that Eq = 
mfa{Hi^) < 0. Define H^^ = H + Vaj as in (4.2) and let Ns{H^,E) he the integrated density of surface states 
of H^ and N{H^, E) the integrated density of states for Hi^. Then 

^.^ ln|lnjV,(g^,^)| _ ln|lniV(ff^,^)| 

E\Eo \n{E - Eq) e\Eo \n{E - Eq) 

where the equality should be interpreted in the following way: if one of the sides exists so does the other one 
and they agree. 

In other words, the Lifshitz exponent for the density of surface states of Hi^ and and the Lifshitz exponent 
for the density of states for H^ agree. 

4.2. Lifshitz tails. In this section we investigate the integrated density of surface states Ns{E) for the 
operator H^ = H + acting on ^^(z^i x'L'^^). We assume throughout that £^0 = 'va.ia{H^^) is (strictly) 
negative and E' < 0. 

By the previous section the investigation of Ns{E) for E near Eq can be reduced to estimates for 
the integrated density of states N{E) of the operator H^^ = H + which acts on ^^(Z'^i). Hence the 
problem of surface Lifshitz tails boils down to ordinary Lifshitz tails in a lower dimensional configuration 
space. However the (free) operator is somewhat more complicated, in fact in Fourier representation it is 
multiplication by 

1 



We remind the reader that K;(7i) = ^^71 for 71 G Z'^i and (<j-'7i)^^g2di is a family of independent random 
variables with a common distribution Pq. 

Throughout this section we assume that supp(Po) is a compact set. Moreover, if we set uj^ = inf(supp(Po)) 
we suppose that Pq{[u!-,u!- + e) > C e^) for some k. 

Theorem 4.2. If h has a unique quadratic minimum then 

^.^ ln|ln(jV,(E))| _ di 

E\Eo ln{E -Eq) 2 ■ 

Proof. The theorem follows from [16, 19] and the considerations above. □ 

For dimensions di = 1 and di = 2 we have the following result: 
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Theorem 4.3. Assume that h is not constant. If di = 1 then 

ln\HN,{E))\ ^ _ ln{n{E)) 
E\Eo ln{E - Eo) E\Eo {E - Eq) 

where n{E) is the integrated density of states for H. 
If d2 = 2, then 

^. ln|ln(7V,(^))| 

lim ^ — = —a 

E\Eo ln{E - Eo) 

where a is defined in (1.8). 

Note that n{E) {E - EqY for some p > 0. See [16, 19] for details. 

To conclude this section we consider some examples that fulfill the assumptions of the previous theorems. 
Let us first assume that H is separable, i.e. that 

h{9i,e2) = hi{ei) + h2{e2). 

This is satisfied for example by the discrete Laplacian where h is equal to /iq given in (1.1). The function 
h has a unique quadratic minimum if and only if both hi and h2 have unique quadratic minima (which we 
may assume to be attained at 9i = 02 = 0). 
We will show in the following that the function 

has a unique quadratic minimum in this case as well. Differentiating the function 

'^^'^ = 1 h,i9i) + h2i92)-Eo'''' 

we obtain: 



{hi{9i) + h2{92) - EoY 



so the (possible) maximum of p is at 9i = 0. 
The second derivative at 0i = is given by: 



l(0) + /!2(«2)-E|,)2 

which obviously gives a negative definite Hessian. 

We remark that no assumptions on /i2 were needed; in fact, the above arguments work for /12 = const as 
well. 

The same reasoning also shows that h is not constant as long as hi is not constant. 
So we have proved: 

Theorem 4.4. Suppose h{9i,92) = hi{9i) + /i2(^2) then 

(1) If hi has a unique quadratic minimum, then 

E\Eo ln{E - Eo) 2 

(2) If di = 1 and h is not constant then 

ln|ln(iV,(i;))| ^ _ ^.^ HnjE)) 
E\Eo ln{E - Eo) E\Eo {E-EoY 

where n{E) is the integrated density of states for H . 

(3) // d2 = 2, then 

^. ln|ln(iV,(S))| 

hm — — 7— — — — = —a 

E\Eo ln{E - Eo) 

where a is defined in (1.8). 
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5. The density of surface states for a constant surfacic potential 



In this section, we prove some useful results on the density of surface states for a constant surface 
potential. In some cases, this density may even be computed explicitly (see e.g. [2]). 

The model we consider is the model introduced in Proposition 1.2 namely Ht = + tl IIq where 
H is chosen as in section 1 and t is a real coupling constant. The proof of all the results we now state is 
based on rank one perturbation theory (see e.g. [24]). The main formula that we will use is the following: 
for z M, one has 

(5.1) (ft - - (H - = - .)-'! nS(ff - 

where I{z) is the operator acting on £'^{Z,'^^) that, in Fourier representation, is the multiplication by the 
function I{6i,z) defined in (1.2). 

Formula (5.1) is easily proved if one makes a partial Fourier transform in the (71,^1) variable of H 
and Ht- If one does so, one obtains a direct integral representation for both H and Ht namely 

H = I H{ei)dei and Ht = [ Ht{ei)d9i 

where H{9i) and Ht{9i) (both acting on differ only by a rank one operator, namely, 

Ht{9i)-H{9i)=ml. 

Formulae (5.1) and (1.2) then follow immediately from the well known resolvent formula for rank one 
perturbations that can be found e.g. in [2 I]. 

Proposition 1.3 follows immediately from Proposition 1.2 and formulae (5.1) and (1.2). Indeed, by 
formula (5.1) and the special form of the operator I{z), z is a point in cr{Ht) \ cr{H) if and only if, for some 
9i, one has 

l + tl{9i,z) = 0. 

If we pick z G M below (recall that = mf{a{H)) = inf (/i(R'^)), we see that z G cr(i?j) if and only if 
tl{9i, z) = —1 for some 9i. As, for z < 0, I{9i, z) is a negative decreasing function of z that tends to when 
z — > —00, we see that this can happen if an only if tl{9i,0) < —1 for some 9i. This is the first statement 
of Proposition 1.1. Indeed, the function 9i 1-^ tl{9i,0) is continuous of T^^i except, possibly, at the points 
where h assumes its minimum, and it takes its minimal value exactly at one of those points. 

As, for the second statement, let I{z) := max I{9i,z) and consider the function / : z 1— > 1 + tl{z). 

This function is clearly continuous and strictly decreasing on ] — 00, 0[ and by assumption, it is negative near 
(as 1 + t/oo < 0) and f{z) — > 1 as z — > —00. So, the function / admits a unique zero that we denote by Eq. 
The analysis given above immediately shows that Eq is the infimum of Ht: as 9i I{9i,z) is continuous 
on T'^i that is compact, for some 9i, one has 1 + t/(^i,£'o) = 0. So that Eq belongs to a{Ht); on the other 
hand, for E < Eq, for any 9i, one has 1 + tI{9i,E) > 1 + tI(E) > 0, hence, E (j{Ht). This completes the 
proof of Proposition 1.1. 

5.1. Asymptotics of the density of surface states. The starting point for this computation is again 
formula (5.1). This enables us to get a very simple formula for the Stieltjes-Hilbert transform of the density 
of surface states n* for the pair {Ht,H). Using the Fourier representation and Parseval's formula, one 
computes 

tT{ni[{Ht - zy^ - {H - zy']Ui) 

Jj-ii l + tli9i,z) Jjd2 h(9i,92)-z h{9i,92)-z ^ 
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sr^ f ~* f e-*T2^2d6'2 f e'"'^^^d92 

^ Ird, l + tl{9i,z)jjd, h{9i,92)-z Jjd, h{9i,92)-z ' 



l + tl{9i,z) V2 {h{9i,92)-z) 
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One then notices that 

7x^1 l + tl{di,z) JYd2 {h[6i,e2) - zY dz Jfd^ 

Here, and in the sequel, log denotes the principal determination of the logarithm. 
This immediately yields that the Stieltjes-Hilbert transform of A^* is given by 

/J_^dNl)=[ \og{i + tiiei,z))dei 

\ - - z I Jjd^ 

where I is defined by (1.2). 

It is well known that one can invert the Stieltjes-Hilbert transform to recover the signed measure dN^s (see 
e.g. the appendix of [2.3]). By the Stieltjes- Perron inversion formula, one has 



(5.2) 



^ 1 /■^// 1 .,.A / 1 



dNt(e) = lim / ( ,dNt ) - ( ,dNt ) de 

= lim — / [ [log{l + tI{9i,e + ie)) -log{l + tI{9i,e-ie))]deide. 
s^o+ 2^7r Jq Jjd^ 



Notice that, for e real. 



Im(l + t/(6'i,e + ie)) = te i 



log(l + tl{9i,e + ie)) - log(l + tl{9i,e - ie)) = log 



'^2 (M^i,e2)-e)2+£2 
hence, this imaginary part keeps a fixed sign. So, for 6i € T"'^ , one has 

l + tl{9i,e + ie) 
I + tl {01,6- ie) 

For e G M, one has \ l + tl{6i, e + i£)\ = \ l + tl{9i, e — ie)\. As moreover the imaginary part of l + tl{9i, e + ie) 
keeps a fixed sign, one has 

I log(l + tl{9i,e + ie)) - log(l + t/(^i, e - ie))\ < 27r. 

As T'^i and [0, E] are compact, one can apply Lebesgue's dominated convergence Theorem to (5.2) and thus 
obtain 

E rE 

t 



(5.3) / dNl{e)= / f{ei,e)deide. 

Jo Jo Jf^i 

where 

(5.4) /(0i, e) = hm J- log f ! t " ^ '1 1 = '^^^^'^ + ^ + ^^)) 

where Arg is the principal determination of the argument of a complex number. Notice here that this 
formula is the analogue of the well-known Birman-Krem formula (see e.g. [1, 25]) for surface perturbations. 

We will now compute the asymptotics of /(^i, e) for e small. First, let us notice that we need only to 
compute these for 9i small, i.e. close to 0. Indeed, we have assumed that h takes its minimum only at 0. 
Therefore, as T*^, is compact, if \6i\ > 5, we know that, for some 6' > 0, for all 62^ one has h{9i,92) > S'. 
Hence, if \9i\ > 6, the function I{Oi, z) is analytic in a neighborhood of 0, so that f{9i, e) = for e sufficiently 
small (independent of ^i). So, we now assume that \6i\ < 5 for some (5 > to be chosen later on. 

We now study I{6i, z) for \z\ small. Pick x a smooth cut-off function in 92, i.e. such that xi^2) = 1 if 
\92\<5^ and x(^2) = if j^al > 25^. Write 

(5.5) m,z)=[ d92+r '-^^^ d92. 

Jjd2 h{ei,92) - z Jjd^ h{ei,92) - z 

For the same reason as above, the second integral in the right hand side term is analytic for \z\ small for all 
9i . We only need to study the integral 

(5.6) J{ei,z)= [ d92. 

Jjd2 h{9i,92) - z 
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Therefore, we use the assumptions that is the unique minimum of h and that it is quadratic non-degenerate. 

This imphes, that for (5 > sufficiently small, for \6i\ < S, the function 62 ^ h{9i,62) has a unique 

minimum, say ^2(^1)) that this minimum is quadratic non-degenerate. Let ^12(^1) be the minimal value, i.e. 

^2(^1) = h{9i,92{0i)). Then, the functions 61 ^2(^1) and 9i 1-^ h{9i) are real analytic in \9i\ < 6. 

All these statements are immediate consequences of the analytic Implicit Function Theorem applied to the 

system of equations V e2^{^ii^'i) = 0. 

So, for 1^1 < one can write 

h{9i,92) = h2{9i) + ((^2 - ^2(^1)), Q2(^i)(^2 - ^2(^1))) + 0{\92 - 92{9i)\^) 
where (52(^1) is the Hessian matrix of /i(^i,^2) at the point ^2(^1)- 

We can now use the analytic Morse Lemma (see e.g. [12]) uniformly in the parameter 9i. That is, for some 
5o > small, there exists i?2(0,(5o) C U (the ball of center and radius 5o in T^^^) ^nd ip{9i) : 92 ^ U ^ 
''/'(^i) ^2) G -62(^2(^1)) 2(5o)) a real analytic diffeomorphism so that, for 9 & U, 

(5.7) h{9i,^p{9i,92)) = h2{9i) + {92, Q2{9i)92). 

Moreover, the Jacobian matrix of ip at ^2(^1) is the identity matrix, and the mapping 9i 1— > ip{9i) is real 
analytic (here, we take the norm in the Banach space of real analytic function in a neighborhood of 0). 
Before we return to the analysis of J, let us describe /i2(') and ^2(0 more precisely. Let Q be the Hessian 
matrix of h at 0. As h has a quadratic non degenerate minimum at 0, Q is definite positive. We can write 
this d X (i-matrix in the form 



(5.8) Q 



Qi R* 
R Q2 



where Qi^2 is the restriction of Q to R'^i''^^ when one decomposes R'^ = M'^i x Both Qi and Q2 

are positive definite; actually, the positive definiteness of Q ensures that the matrices Qi — R*Q2^R and 
Q2 — RQi^R* are positive definite. Using the Taylor expansion of h near 0, one computes 

^2(^1) = -Q2'R0i + 0{\9i\^), Q2{9i) = Q2 + 0{\9i\), 

(5.9) h2i9i) = {[Qi - R*Q^^R]9i, 9i) + 0{\9i\'). 
Let us also note here that 

(5.10) Det Q = Det Qi • Det (Q2 - R*Qi^R) = Det Q2 ■ Det (Qi - RQ^^R*) 
We now return to J. Performing the change of variables 9 ip{9) in J{9i,z), we get 

(5.11) Ji9i,z)= [ r/lnl^yu ,n ^ d^2wherex(^i,02):=x(V'(^i,^2))Det(Ve,m,02)). 

Choosing 6 sufficiently small with respect to 5^ (defining %), we see that x(''/'(^i) ^2)) = 1 for all |^i| < 6 
and 1^2! < Hence, the function x(^i)^2) is real analytic in a neighborhood of (0,0). 

To compute the integral in the right hand side of (5.11), we change to polar coordinates (recall that x 
is supported near 0) to obtain 

I-+00 ^(n \ d2-l 

(5.12) J{9,,z) = Det (Q2(^i))-^/' / 2 1 L , dr 

Jo H + ^2(6*1) - z 

where 



(5.13) x{Oi,r):=-—^ m,rm- 



-1 



The factor (27r)~'^2 in the last integral comes from the fact that d92 denotes the normalized Haar measure 
on ¥"^2^ i_g_ Lebesgue measure divided by (27r)'^2. Note again that (^i,r) 1-^ x(^i)''') is real analytic in a 
neighborhood of 0, and 

X(^i,0) = (^Det (V,,V(^i,^2(^i))) • Vol(S^2-i). 
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Moreover, as f^d2 ^^d^ = if A; is multi- index of odd length, we known that the Taylor expansion of 
contains only even powers of r, i.e. there exists a function x(0i,r) analytic in a neighborhood of (0, 0) such 
that xiOi,r) = x{ei,r'^). 
We now use the 

Lemma 5.1. Let x he a smooth compactly supported function such that x be real analytic is a neighborhood 
ofO. Define the integral Jj^{z) to be 

M^) = I ^:2^dr. 

Then, one has 

(5.14) J^{z) = S{z) ■ H{z) + G{z) 
where 

(1) G and H are real analytic in a neighborhood of 0; 

(2) they satisfy H{0) = x(0) and G(0) > if x(0) > and x > 0; 

(3) the function S is defined by 

• if n is even, then S{z) = — ■ [—l)2z 2 . logz; 

^ / s Hni "-1 1 

• if n is odd, then S{z) = — • (—1) 2^2 ——, 

2 

Here, \fz and logz denote respectively the principal determination of the square root and of the 
logarithm. 

The proof of this result is elementary; after a cut-off near zero, one expands x 'va. a, Taylor series near 0, and 
computes the resulting integrals term by term essentially explicitly (see [15] for more details). 
Putting (5.5), (5.6), (5.11), (5.12) and (5.14) together, we obtain that 

(5.15) i{ei,z) = s{h2{9i) - z) ■ H{9i, h2{9i) - z) + G{eiM{ei) - z) 

where 

• S* is described in point (3) of Lemma 5.1; 

• (^1, z) I— > H{6i, z) and (^i, z) 1— > G{Oi, z) are real analytic in 9i and z in a neighborhood of 0; 

• one has 

m,0) = ^^Det (Q2(0i))"'/'-Det (V,,V'(^i, e2))Vol(S'^^-i) 
and G(0, 0) is positive. 

The last point here is obtained combining point (2) of Lemma 5.1, (5.12) and (5.13), and using the decom- 
position (5.5). 

The first immediate consequence of (5.15) is that, if e € M and /i2(^i) > e, then 

I{ei, e + ie)- I{9i,e - ie) ^ when e ^ 0+. 
This implies that, if /i2(^i) > e, one has 

/(^i,e) =0. 

Assume now that h2{0i) < e. As < h2{6i), —e < /i2(^i) — e < 0. We now need to distinguish different 
cases according to the dimension d2- Consider the case 

• d2 = 1- by (5.15), as H and G are analytic, one has 

hm 1(01, e + ie) = ' =g(gi, /t2(gi) - e) + G(0i, /i2(ei) - e). 

e^0+ I \/\h2\0\) - e| 

Using again the fact that H and G are analytic and that H{6i,0) does not vanish for 61 small, we 

get 

l + tl{9i,e + i£) _ ^ , .2G(0,0)y|/i2(6'i) -e| 



lim — '-^ — ■ = -1 + i ^ ' \ ' + o J\h2(ei 

e^o+ l + tl{ei,e-ie) tH{0,0) iviAi 
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As G(0, 0), H{0,0) and t are also positive, one finally obtains 



/(^i,e) = i[l + o(V|/i2(ei)-e| 



1{/J2(9l)<e}- 



d2 = 2: in this case, one computes 



lim I{ei, e + ie) = ^(1 log |/i2(0i) - e|| + iir)H{ei, /izl^i) - e) + ^(^1, /i2(^i) - e). 

e^0+ 2 

Using again the fact that H and G are analytic, we get 

i + ^ A J';; ) . d + oioog - ei)-i). 

e^o+ 1 + t/(6'i, e - ze) y I log 1^-2(^1) - e|| / 
So that finally, one has 

= |log|/.2(^gi)-e|| ^^ + '^^^^^^ " "'^"'^^ ■ 

• £^2 > 3: in this case, one has to distinguish two cases whether 1 + tl{0, 0) = or not, as well as the 
case of even and odd dimensions. 
Let us first assume: 

— that 1 + 1/(0, 0) > 0: as H and G are analytic, one has 



lim 1(01, e + ie)= { lim 5(/i2(^i) - e - ie) ] ■ HiOi, h2(9i) - e) + G(0i, /laf^i) - e). 

As G is analytic and as /^(O) = 0, one has G{9i,h2{0i)) = I{9i,0). So, for 9i small, we know 
that 1 + tG{9i, 0) 7^ 0. Here, we used the continuity of G and the fact that ^2(^1) is of size j^ip 
hence small. This gives 

l + tl{9^,e + ie) _ t-s{h2{9i)-e)-H{9i,0) 
/^^+l + tl{9i,e-ie) ^ l + t-G{9i,0) ' + 

where 

s{x) = lim [S{x - ie) - S{x + ie)], = O((/i2(0i) - e) • \S{h2{9i) - e)|, (/i2(^i) - e)). 
So that finally, for e small, one has 

f^ifi^ f(ft ^ t.sih2{9,)-e)-H{9uO) , ^ 
(5-16) fi9i, e) = i + tG{9i 0) Mh2{e^)<e}{l + 

where 

1 

(5.17) s{x) = -\x\ 2 

and R is given above. 

Prom these asymptotics and from (5.3), integrating / in (5.3), using (5.9) and (5.10), one gets that 
• if (i9 = 1: 



(5.18) 

• if d2 = 2: 



•^0 di{di + 2)(2^)'^i -y/Det (Qi - RQ^^R*) 



^,^.(,) _ 2Vol(S^^-) E^^^^/^ 



di(di + 2)(27r)<^i -^Det (Qi - RQ^^R*) 
• if d2 > 3 and 1 + t/(0, 0) > 0: 

(5.19) frn'M ~ °(^i.''2)yo'(s;-;^(s^-') 

^ ^ Jo E-.0+ (i(27r)'^VDetQ l + t/(0,0) ^ ^ 
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where S'^^'^ ^ are respectively the di^2 — 1 dimensional unit spheres, and s is given by (5.17). Here, c{di,d2) 
is the integral 

(5.20) c(di,d2)= rr'^i-i(l-r2)('^2-2)/2^^_ 



5.1.1. The borderline case. Though it will not find direct applications in this paper, let us now turn to the 
case when c?2 > 3 and 1 + tl{0, 0) = 0. Notice that this assumption implies t < 0. When 1 + tl{0, 0) = 0, 
one has to take a closer look at the vanishing of 1 + tG{6i, 0) when 9i 0. We will now assume that 

(H): I (01,0) has a local maximum for 9i = 0. 

Remark 5.1. Notice that this assumption was also necessary for fluctuating edges. Actually, in that 
setting, we even required that the maximum be non-degenerate if di > 3. This seems quite natural as the 
case 1 + tl{0, 0) = is exactly the border line between the fluctuating edges and stable edges. 

Let us recall that as above, we need to compute the asymptotic when e — > 0"^ of the integral 

/ f{ei,e)dei= [ f{ei,e)dei 

JT'^i J{h2(ei)<e} 

where / is defined by (5.4). Using (5.9) we can find a analytic change of variable 6i ip{Oi) such that 
/i2(^~n^i)) = {QiSi,Oi) =: 92(6*1) where Qi = Qi- R*Q2^R (the matrices (5i,2 and R are defined in (5.8)) 
and ^(^i) = ^1 + 0(1^1 p). So, we want to study 

f{iP{9i),e)\DetVeM(^i)\dei 

{q2{0i)<e} 

Let us perform one more change of variable in the integral above, namely 0i y/eOi; hence, we need to 
study 

[ f{^P{V^6i),e)\DetVeMV^^i)\d0i 
J{Q2{ei)<i} 

Notice that, for e small, on {{Qi0i,Oi) < 1}, one has 

|DetV0,V(V^^i)l = l + OiV^). 
We now study /(V'(\/e^i), e) for e small and {^2(^1) < 1}- 

Using the analyticity of G and H, for e > 0, we start with rewriting (5.15) in the following way 

(5.21) l + tI{'^iV^ei),e + ie) = 1 + tG(?/'( V^^i), 0) + te5,G(0, 0)((?2(^i) - 1) + 

+ tSie ■ ((72(^1) - 1) - ie))H{^P{V^9i),0) + 0{e + + |e • S{e)\). 

Let us now distinguish between the different dimensions, i.e. between the cases d2 = 3, d2 = 4: and ^2 > 5. 
Substituting the asymptotics for S given in Lemma 5.1 and using the analyticity of G and H, one obtains 
the following: 

• If d2 = 3: define F±(6i,e) = lim 1 + tlNji^/eOi), e ± ie). For ^2(^1) < 1, one has 
F±(0i, e) = (Tit(2^)-''^|1 - g2(^i)|'/'Det (Q2)~'/' + t ■ g{ei) + o{V^)^ 

where 

(5.22) t-g{ei)= lim ^[l + tG{^P{V^ei),0)]. 

e^0+ -y/e 

This gives, for 92(^1) < 1, 

/(V(Viei),e) ~ iArg(-i(27r)-'^^Det(g2)-'/'|l- '72(^1)1'/' + 5(^1)). 
e^o+ vr 
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We notice that this last argument is non positive. As a result we obtain that 

[ Avg{-i\i-ej\^/' + ~g{e,))dei 

(5.23) / dNl{e)de ~ = • 

•^0 di{di + 2)7r(2^)'^i -^Det (Qi - RQ^^R*) 

where 

= (2^)'^VDet (Q2)5((Qi - /?Q2'^*)"'/'^i) 

and g is defined by (5.22). 

Remark 5.2. In some cases, g and (7 are identically vanishing. This happens for example if /i is a 
"separate variable" function, i.e. if h{6i,62) = /ii(^i) + /i2(^2)- Indeed, in this case, /i2(^i) = ^i(^i) 
and I{6i, h2{6i)) = 1(0,0), hence, G does not depend on 9i, i.e. G{6i,z) = G{z). 
When g vanishes identically, formula (5.23) becomes (5.18) except for the sign which changes to — . 

The integral / Arg(— — + g{6i))d9i is negative. Hence, comparing (5.23) to (5.19), 

J\ei\<i 

we see that, asymptotically when E 0^, Jq dN^{e)de is larger when 1 + t/(0,0) = than when 
1 + tl{0, 0) > 0. This is explained by the fact that, when 1 + tl{0, 0) = 0, a zero energy resonance 
(or eigenvalue if ^2 > 5) is created. This resonance (eigenvalue) carries more weight. Of course, the 
same phenomenon happens for the spectral shift function. 
To conclude the case d2 = 3, let us notice that we did not use assumption (H). 

• If d2 = 4: let us start with computing dej^G{0,0). Therefore, we use G{6i,0) = I{6i, h2{0i)) and 
compute 



90i/l(O,02) 



|ei=o 



:dei = 



Ij,, (/l(O,02))2 

as is a local maximum of I{9i,0). This computation immediately gives that 1 + tG{'^p{^/e9l), 0) 
0(e). Hence, equation (5.21) gives 

F±{9ue) = te{q2{9i) - l)(log e + log 1^2(^1) - 1| + h{e)) • 1 + 



log e + log \q2{9i) - 1| + h{e) 
where h{e) is bounded and does not depend on the sign it. This gives, for (72(^1) < 1, 

/(V^(V;0i),e) ~ 

e^o+ I log e I 

Integrating over 9i and e, we obtain 

^^^^ 2Vol(S'^i-i) 



di {di + 2) (271)^^1 JDet {Qi-RQ^^R*) ' ^1 



dNlie) 

If 1 + f/(0, 0) = and ^2 > 5: we now compute dg^Q{0,0). Therefore, we continue the computation 
done above to obtain 



|ei=o 
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d'g^{h{9i,92)-h2{9i)) \ f f [deAhi0i,92)-h2i9i))] 



^^•^^^ (M^i,^2)-/i2(ei))2 "^^Vl^^^/H^. {H9i,92)-h2{9^)r '^^Vi.,=o 

''^^'^'^''^d9,^2f ^W^d9,J[ ^^d9,]Q2 



T<^2 (/i(0,^2))2 M {h{0,92)f \Jj,, {h{0,92)) 
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where Q2 is defined in (5.8). On the other hand, one has 

a,G(0,0) =-a,/(0,z)|,=o = -^ where J := / ^ d9i. 

Plugging this and (5.24) into (5.21), we obtain 

1 + t/(V'(^/^^i), e ± ie) = -teJ + o(e) + tS{e ■ (^2(^1) - 1) T ie){H{0, 0) + o(l)). 
where o(e) does not depend of ±. This gives, for 92(^1) < 1, 

.(e • (92(^1) -1)) 



*o+ J 

Integrating over 61 and e, we obtain 
where c{di,d2) is defined in (5.20). 

6. Appendix 

Pick E < —d. We now prove that, for h taken as in Remark 1.1, the function h defined in (1.7) is not 
constant. For the purpose of this argument, we write 9i = {6^, ■ ■ ■ ,6'^^). To check that h is not constant, 
by (1.2) and (1.7), it suffices to check that the function 61 ^ I{6i,E) is not constant, hence, that the 
function 9^ Ji(^^) defined by 

(6.1) j{e') = —L_^ [ ii^d\e\... y\E)de^...de^- = -\^^ [ i de' 

is not constant. We used the notation 6 = (^1, 62) = {0^, 9'). 

Recall from Remark 1.1 that h{6) = /io(G' • 9) where G' G GSL(Z) and /iq is defined in (1.1). So, the n-th 
Fourier coefficient of J is given by 

J = ^ [ [ .p^.pl = _}_ f ,n 

" (2vr)'^ Jlo,2.] y[o,2.]^-i h{9\ 9') - E (27r)^ jy,^2.Y ^o(G' ■ 9) - E 

-d9 = — / — -d9 



{27tY Jio,2n]^ ho{9) - E {27ry -ho{9) - E 

where {G'~^ ■ 9Y denotes the first coordinate of the vector G'"^ ■ 9, and 9-,^, the vector (vr, . . . , vr) in W^. So 
to prove that J„ does not vanish for any n which implies that J is not constant, it suffices to prove that the 
Fourier coefficients of {hQ{9) — E)^^ do not vanish. This is a consequence of the Neuman expansion 

1 -I fho{9y'' 



-ho(9) -E E ^ \ -E 

Indeed, the n-th Fourier coefficient in each of the terms of order k larger than n in this series is positive 
: it is easily seen as —E > and the multiplication operator (/iq)" is unitarily equivalent through Fourier 
transformation to (— ^A)"; so the Fourier coefficients of (/iq)" are the entries of the zeroth row of the matrix 
(— ^A)" and, the n first super- and sub-diagonals of this convolution matrix are positive. 
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